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Abstract 

We study three and four jet production in hadronic collisions at next-to-leading order accuracy in massless QCD. We 
cross check results previously obtained by the BlackHat collaboration for the LHC with a centre-of-mass energy 
of yfs — 7 TeV and present new results for the LHC operating at 8 TeV. We find large negative NLO corrections 
reducing the leading-order cross sections by about 40-50%. Furthermore we observe an important reduction of the 
scale uncertainty. In addition to the cross sections we also present results for differential distributions. The dynamical 
renormalization/factorization scale used in the calculation leads to a remarkably stable K-factor. The results presented 
here were obtained with the NJet package [ 1 1, a publicly available library for the evaluation of one-loop amplitudes 
in massless QCD. 
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1. Introduction 

Multi-jet production in hadronic collisions via the 
strong interaction represents an important testing 
ground for quantum chromodynamics (QCD). Jet pro- 
duction in massless QCD may provide valuable in- 
formation to constrain the QCD coupling constant oc s 
and/or the parton distribution functions (PDFs). Fur- 
thermore owing to their significant production rates 
multi-jet processes may contribute as background pro- 
cesses to a large variety of interesting signal reactions 
including processes relevant for new physics searches. 
Independent of whether multi-jet production in QCD is 
studied as signal or background process precise theoret- 
ical predictions are mandatory. 

In the past substantial progress has been made in 
the numerical evaluation of leading order (LO) ma- 
trix elements for large multiplicities. Using publicly 
available tools like Alpgen [2|, Sherpa+Comix J3] O, 
MadGraph/MadEvent multiplicities with up to 12 
final state jets can be calculated. However leading- 
order cross section predictions suffer from a substantial 
(residual) dependence on the unphysical renormaliza- 



tion scale /j r . This is in particular true for high multi- 
plicities due to the high power in oc 4 occurring in the 
perturbative expansion. For a process starting with a" 
the scale dependent terms are proportional to 

n$ a s O (1) 

where Gq denotes the leading-order cross section and 
Po is the leading coefficient of the QCD beta function. 
Precise theoretical predictions thus require at least next- 
to-leading order (NLO) corrections in the coupling oc s 
to reduce the scale dependence. While the NLO cor- 
rections to two jet production were calculated already 
in 1992 J8l|9l and the corrections to three jet produc- 
tion a decade later ifTOlfTT l (results for the gluon channel 
were obtained already in Refs. lfl2l[T3ll ) further progress 
has been hindered by both the increasing complexity of 
the one-loop amplitudes and the large number of tree- 
level processes entering the infra-red counter-terms in 
the Catani-Seymour subtraction formalism [14|. The 
latter problem could be solved by an automated gen- 
eration of subtraction terms based on the aforemen- 
tioned tree-level technology flSlEOl . While techni- 
cally solved the available computer resources still im- 
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pose a severe restriction on the multiplicity which can 
be handled in practice. As far as the virtual corrections 
are concerned just the mere number of Feynman dia- 
grams in the conventional Feynman diagram based ap- 
proach e.g. ~ 15, 000 for gg — > 4g at one-loop (neglect- 
ing self-energy type corrections) shows that alternative 
techniques are required. In recent years new methods 
for the evaluation of one-loop amplitudes have relieved 
the long standing bottleneck in producing virtual correc- 
tions for processes with high parton multiplicities. The 
on-shell unitarity method 12T1 l22l and the subsequently 
developed generalized unitarity method [23 -28 1 offer 
an interesting alternative to conventional Feynman dia- 
gram based methods. In particular the unitarity method 
allows to use leading order amplitudes as basic build- 
ing blocks for the one-loop calculation. It is thus pos- 
sible to make use of the efficient tree-level machinery 
mentioned above. In parallel with the development of 
the unitarity method also the conventional approach has 
been continuously improved in the last couple of years 
and a mixture of methods are currently exploited |29l . 
We also mention that recently other directions have been 
investigated which offer the potential to become pow- 
erful alternatives to those above ll30l . While improve- 
ments are still ongoing, NLO technology is already now 
automated and flexible enough to cover a wide range 
of processes for four final state particles PTH491 and a 
handful of 2 ->> 5 processes (to] I50H521 . A number 
of automated approaches to the computation of virtual 
amplitudes have appeared as public codes ll53Tl56l . 

Despite the aforementioned progress only very re- 
cently first results for four jet production in NLO ac- 
curacy became available. In Ref. Il57l the Black- 
Hat collaboration published results for three and four 
jet production in NLO accuracy for the LHC running 
at 7 TeV centre-of-mass energy. Given the complex- 
ity of the calculation it is very important to perform an 
independent cross check. Furthermore we present re- 
sults for the LHC running at 8 TeV and discuss vari- 
ous differential distributions. For the evaluation of the 
virtual matrix elements we use the publicly available 
package NGLUON (53]. NGLUON uses on-shell meth- 
ods to evaluate numerically one-loop primitive ampli- 
tudes in pure gauge theory. Recently we have extended 
NGLUON to account for amplitudes involving massless 
quarks CD ED ED- We note that in Refs. CD ED ED no 
approximation in the colour is performed. The real cor- 
rections and the cancellation of soft- and collinear diver- 
gences are obtained within the Catani-Seymour dipole 
subtraction method |14|. For this contribution we use 
the Sherpa Monte-Carlo event generator ifTBI . We will 
briefly comment on some technical aspects in the next 



section. In section 3 we will present the comparison 
with the results published by BlackHat ||57). In addi- 
tion we show new results for 8 TeV for three and four jet 
production. Our conclusions are presented in section 4. 

2. Outline of the calculation 

In what follows we briefly outline the calculation. We 
work in massless QCD with 5 flavours. In particular 
we include the bottom quark in the initial state. The 
top-quark is assumed to be integrated out. Through the 
matching of the coupling constants between the five and 
the six flavour theories some of the corrections due to 
top-quarks are retained. We expect the neglected cor- 
rections to be small. For three jet production the con- 
tributing processes are given by all possible crossings 
of the following transitions: 

O^ggggg, Q->qqggg, O^qqq'qg, (2) 

where q and q' denote a generic quark. For the four 
jet production the corresponding processes are derived 
from 

-> gggggg, (3) 

and 

O^qqgggg, Q^qqq'q'gg, 0^ qqqqq'q" '. (4) 

In principle we need to distinguish the like-flavour cases 
from the case that all quark flavours are different. How- 
ever technically the like-flavour processes can be ob- 
tained from the latter by an appropriate (anti) sym- 
metrization with respect to the quark momenta. For ex- 
ample the amplitude A q qqqg for the transition — s- qqqqg 
is obtained through 

^gfg3j(l)2,3,4,5) = 

^i?ot'^(1j2,3,4,5) — A^Yg(l,4,3,2,5), (5) 

where we used the short hand notation to abbreviate 
with i the momentum kj and helicity A., of the z'-th par- 
ton (i = (ki,Xi)). The expansion of the n-jet differential 
cross section in the coupling a s reads: 

da„ = da\;° + 8do?, LO + 0(a"+ 2 ) (6) 

(c/aj; ~ a",5c/a^ L0 — a" +1 ) with the leading order 
cross section given by 

d ^n°= E dx i dx j F j/H 2 (xj,/Jf)Fi/H l (xi,Vf) 

£{q,q,g} 

x dol (i(xiPx ) + j{ X jP 2 ) -> n pail.) . (7) 
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Pi, Pi denote the momenta of the incoming hadrons. 
Xj are the momentum fractions carried by the initial 
state partons i with respect to the incoming hadrons 
H\,H%. The momentum sum of the incoming partons 
is thus given by P = XiPi +XjP2- As usual we ne- 
glect the masses of the incoming hadrons. F i / H (xi,^if) 
are the parton distribution functions which roughly 
speaking describe the probability to find a parton i 
inside the hadron H with a momentum fraction be- 
tween Xi and Xi + dxj. The factorization scale is de- 
noted by /jf. The partonic cross section for the reac- 
tion (ij — > n-jets) in Born approximation is given by 
dafj (i(xiPi) + j(xjPi) —> n-jets). In terms of the lead- 
ing order squared matrix elements |iW„| 2 the explicit ex- 
pression for dafj reads: 

B _ 1 A d ^ k t 

x (27t) 4 5 (P - £ * m ) | <M n (ij -> n part.) | \ (8) 

m= 1 

where kj are the four momenta of the outgoing partons 
and S — 2xfXj(Pi -Pi) = (L"=i^i) 2 * s tne partonic cen- 
tre of mass energy squared. The jet algorithm is en- 
coded through 0,,-jet- ©n jet is a function of the final 
state parton momenta kj. It is one if the corresponding 
n parton momenta represent an n-jet event and zero oth- 
erwise. As mentioned in the introduction different pub- 
licly available tools exist to evaluate the matrix elements 
M(ij — > n part.). We used Amegic++ |60| within the 
Sherpa framework to achieve this task. Cross checks on 
the matrix elements were performed using Comix J4|. 
To perform the numerical integration the Sherpa Monte 
Carlo event generator lfl5l is used. 

At NLO accuracy we need to consider the one-loop 
corrections da^ and the real corrections daf +l due to 
an additional parton in the final state. Both contribu- 
tions da^ and daf +i individually contain collinear and 
soft singularities. To obtain a finite result the two contri- 
butions must be combined and the initial state singulari- 
ties must be factorized. A convenient method to perform 
the cancellation of the soft and collinear divergences is 
the Catani-Seymour subtraction method fl4l . The ba- 
sic idea is to introduce local counter-terms which render 
the integration of the real corrections finite and are easy 
enough so that an analytic integration over the soft and 
collinear regions of phase space is possible. Schemati- 



cally the total cross section reads: 
5o NLO = f (d^ + J da s n+1 ) + J da**- 

n 1 n 

+ J {da% +1 -da s n+l ), (9) 

n+l 

where da^ +l denotes the local counter-term and da^ ac - 
is due to the factorization of initial state singularities. 
It is thus convenient to split the NLO corrections into 
three contributions: 

8da^ L0 = do^ + dd\ + da^ , (10) 

the finite part of the virtual corrections do^, the finite 
part of the integrated subtraction terms together with 
the contribution from the factorization d<5 l n and the real 
corrections combined with the subtraction terms rfcr^fj . 
As in the case of the LO cross sections we use Sherpa in 
combination with Amegic++ to calculate da], and daf s . 
Again cross checks on the tree-level amplitudes were 
performed using Comix [0). 

The one-loop matrix elements required for the virtual 
corrections da^ are evaluated using an on-shell gen- 
eralized unitarity set-up (see Ref. |61 1 for a recent re- 
view) for multi-parton primitive amplitudes [53 , 58l l59l . 
In Ref. [53 1 only colour-ordered pure gluon amplitudes 
were considered. To account for massless quarks ap- 
pearing in the loop or as external partons we have ex- 
tended the NGluon package to allow the computation 
of primitive amplitudes involving quarks. The primitive 
amplitudes correspond to the colour-ordered amplitudes 
in the pure gluon case. Details on the extension of the 
NGluon package are given in Refs. |H||58]I22]|- The 
entire library together with additional code to perform 
the colour algebra — on which we comment below — 
is publicly available as NjElQpack age. A detailed de- 
scription how to use the library can be found in Ref. QJ. 

For a given process the primitive one-loop ampli- 
tudes provide all required information to reconstruct the 
full amplitude including the complete colour informa- 
tion. In particular the partial amplitudes appearing in 
the colour decomposition of the full amplitude can be 
obtained as linear combinations of primitive amplitudes. 
Since in our calculation no approximation in the sum 
over colour is performed, we need to express all partial 
amplitudes in terms of the primitive amplitudes. This is 
in general a non-trivial task. In Ref. ll6D an algorithm 
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to establish the relation between the partial amplitudes 
and the primitive amplitudes has been presented. The 
method uses the colour decomposition of the full am- 
plitude in terms of colour stripped Feynman diagrams. 
This decomposition is obtained by separating for each 
Feynman diagram the Lorentz structure from the colour 
structure. On the other hand the primitive amplitudes 
can also be expressed as linear combinations of colour 
stripped Feynman diagrams. From the comparison of 
the two representations the relation between the partial 
and the primitive amplitudes can be extracted. For ad- 
ditional details on the method we refer to Ref. ll6D . We 
have applied this method to produce results for up to 
seven point amplitudes ll59l . The explicit formulae have 
also been implemented in the NJet library. We note that 
the relations between the primitive amplitudes and the 
partial amplitudes for up to seven partons are also given 
in Ref. Il62l . Since we apply additional symmetries we 
slightly differ in the number of independent primitive 
amplitudes required in the numerical evaluation. To 
check the correctness of our implementation of the vir- 
tual corrections we have compared our results for indi- 
vidual phase points as far as possible with GoSam ll56l 
and Helac-NLO ll55l . Furthermore we also checked 
the benchmark points provided by BlackHat ||57ll . 
We obtained at least an agreement of eight digits. This 
is largely sufficient for all practical applications. 

To perform the phase space integration we interfaced 
the NJet library to the Sherpa event generator using the 
Binoth Les Houches Accord ll63l . As a technical remark 
we add that we used the HepMC file format f64l for 
weighted events. 

3. Cross sections for multi-Jet production at the 
LHC 

3.1. Numerical setup and checks 

Before discussing the results for three and four jet 
production at the LHC in massless QCD we briefly de- 
scribe the numerical setup adopted in the calculation. 
As shown in Eq. <JVj we need to specify the parton distri- 
bution functions. We use the MSTW2008 PDF set |§3. 
The MSTW2008 PDF set also provides a set of error 
PDFs which can be used to assess the uncertainties due 
to our incomplete knowledge of the PDFs. We have not 
performed any detailed analysis of PDF uncertainties 
or comparison between alternative fits Refs. ||66H68l , 
which we leave for a future publication. In particular 
it would be interesting to see how the MSTW2008 PDF 
set compares with the ABM12 set [68], which comes 
with a slightly different value for a s and differs signifi- 
cantly in the gluon distribution at large x. 



For consistency we use the a s values as provided by 
the PDF set. We note that the MSTW2008 PDF set 
contains so-called leading-order, next-to-leading order 
and next-to-next-to-leading order PDFs. They were ob- 
tained by using different orders in the evolution and 
different orders for the hard scattering coefficients in 
the fits to data. Since U s is extracted together with 
the PDFs different orders come in general with differ- 
ent oCi values. In particular for the LO set we have 
a^°(jj r = m z ) = 0.13939 while for the NLO set the 
corresponding value is a^ LO (p r = mz) = 0.12018. The 
large value of oc^° reflects the fact that LO PDF fits in 
general don't work very well since sizeable NLO cor- 
rections are not taken into account. This in turn leads to 
the large value for oc 4 which partially compensates the 
missing corrections in the hard scattering coefficients. 
Despite this obvious tension with the world average for 
a,, we follow in our default setup the standard procedure 
to use LO PDFs together with ajr° to evaluate LO cross 
sections. For the NLO cross sections we use the NLO 
PDFs together with a^ LO everywhere. When we dis- 
cuss the size of the NLO corrections we will come back 
to this point. The PDFs and 0C S dependent on the un- 
physical scales /j/ and p. r . For the results presented here 
we set p.f = fi r = /j. For distributions where the typi- 
cal energy scale Q changes significantly like for exam- 
ple the transverse momentum distribution of the highest 
energetic jet a fixed value of /j may lead to a poor be- 
haviour of the perturbative expansion due to the appear- 
ance of possible large logarithms of the form ln(/j/Q). 
In such cases using a phase space dependent a s is usu- 
ally a better choice since some of the logarithms may 
be resummed through the evolution of a s . As a conse- 
quence we use a dynamical scale /j based on the sum of 
the total transverse momentum of the final state partons 

^parton 

h t = E p p r a ; ton - (id 

In particular we set /j = Hj /2. To estimate the effect of 
uncalculated higher orders we consider the scale varia- 
tion H T /4 </j <H T - 

For the jet algorithm appearing in Eq. ([8]) through 
©n-jet we use the anti-kt algorithm |69| as implemented 
in FASTJET [70|. The jet-radius parameter/? is set to 

R = 0A (12) 

following the value adopted by the ATLAS collabora- 
tion. Events were generated using identical cuts to that 
of the multi-jet measurements from ATLAS ll7TI and the 
recent study by the BlackHat collaboration |57]]. In 
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particular the transverse momentum, pj, of the first jet 
is required to be larger than 80 GeV with subsequent jets 
required to have at least pj > 60 GeV. Rapidity cuts of 
|t|| < 2.8 were also taken. 

To check our setup we reproduced the results quoted 
in Ref. Il57l . For the three and four jet cross sections at 
a centre-of-mass energy of 7 TeV we find: 



,7TeV-LO 



,7TeV-NLO 



: 93.40(0.03)^4^ 



= 53.74(0.16)^6^ 



20.72 " u ' 



and 



,7TeV-LO 



= 9.97(0.02) nb, 



„7TeV-NLO 
°4 



= 5.56(0.17) nb, 
where the NLO cross section 0^ LO is defined as 

8a? L0 . 



_NLO _LO 

o„ = a„ 



(13) 
(14) 

(15) 
(16) 

(17) 



The number in brackets show the Monte-Carlo errors 
of the numerical integration and the sub- and super- 
scripts refer to the minimum and maximum values ob- 
tained through scale variations estimated at p. = Hj /4 
and p — Hj. In case of the four jet rate we have only 
compared the central scale to reduce the computational 
effort. Comparing the results given in Eq. (14i and 
Eq. ( 16 1 we find perfect agreement with Ref. 0571 . We 



believe that together with the check of the matrix ele- 
ments for individual phase space points the successful 
comparison with the Ref. {571 is highly non-trivial giv- 
ing us confidence that the results presented in this paper 
are correct. 

3.2. Three Jet Production 

We now present results for three jet production at 8 
TeV centre-of-mass energy. Since three jet production 
has been studied in some detail before ifTOlfTTI we limit 
our discussion to a few basic quantities. Using the afore- 
mentioned setup we find for the total three jet cross sec- 
tion: 



,8TeV-LO 



„8TeV-NLO 
°3 



= 126.65(0.05)+^ nb, 
= 72.57(0.16) 



+2.71 
-28.08 



nb, 



(18) 
(19) 



again the numbers in parentheses quote the numerical 
uncertainty due to the Monte Carlo integration while the 
sub- and super-scripts show the effect of the scale varia- 
tion. We note a significant reduction of the scale uncer- 
tainty similar to what has been observed for 7 TeV. Fur- 
thermore the NLO corrections give a sizeable change 
of the cross section prediction: The NLO results are re- 
duced by about 40% with respect to the LO cross section 
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Figure I: pr distribution for the leading jet at the LHC with a centre- 
of-mass energy of 8 TeV. The upper plots show leading order (blue) 
and next-to-leading order (red) results for the central scale p = Hj /2. 
The bands show the LO and NLO scale variations respectively. In the 
lower plot we show the ratio of LO and NLO together with the LO 
scale variations (blue band) and NLO scale variations (red band). 



again in perfect agreement to what has been observed 
for 7 TeV in Ref. l57l . Compared to a collider energy 
of 7 TeV the cross sections are about 35% larger as a 
consequence of the increased parton fluxes at 8 TeV. As 
a technical side remark we point out that the numerical 
integration is very well under control: The numerical 
uncertainty for the central value amounts to about 2 per 
mille. In Fig. [T] we show results for the pj distribution 
of the leading jet. As for the total three jet corrections 
we observe a sizeable reduction of the NLO prediction 
with respect to the leading order result. The scale un- 
certainty is reduced to about 25% of the LO scale un- 
certainty. With the exception of the low pj region we 
observe a rather constant K-factor. The dynamical scale 
chosen for the renormalization and factorization scale 
indeed avoids large logarithmic corrections at large pj. 
For small pj we expect that soft gluon corrections be- 
come important. A reliable prediction in that region 
would thus require to go beyond a fixed order calcu- 
lationj^We have also studied the pj distribution for the 



2 Strictly speaking already the dynamical scale setting procedure 
is beyond a fixed order calculation. A certain class of possible large 
logarithms are resummed via the renormalization group. In case of 
soft logarithms in addition soft gluon resummation would be required. 



5 



T3 



b 

T3 




2.0 
1.5 
1.0 
0.5 
0.0 



NJet + Sherpa 

pp^3jet at 8 TeV 



-10 1 
Leading jet rapidity 



Figure 2: Rapidity distribution in pp — > 3 jets for the leading jet at the 
LHC with a centre-of-mass energy of 8 TeV. 



second and the third jet in the pj ordering. The results 
are very similar to what is shown in Fig. [T] Again we 
find a rather flat K-factor. The main difference to the 
pj distribution of the leading jet is that the corrections 
are slightly larger. The results are shown in the appendix 
(Fig. A. 1 1 and Fig. A. 12 1. In Fig.[2]we show the rapidity 
distribution of the leading jet. As in the pj distribution 
we observe a sizeable correction at NLO together with 
an important reduction of the scale uncertainty. Again 
the K-factor is rather flat. We have also analysed the 
rapidity distribution of the second and third jet. The re- 
sults look very similar to Fig. [2] The explicit plots are 
given in the appendix. 



3.3. Four Jet Production 

In this subsection we present new results for four jet 
production at NLO accuracy in QCD. Performing a sim- 
ilar analysis as in the three jet case we find for the four 
jet cross section at 8 TeV: 



,8TeV-LO 



,8TeV-NLO 



= 14.36(0.01)+^° 6 38 nb 



= 8.15(0.09)+^ 4 nb 



(20) 
(21) 



Similar to the findings for 7 TeV the size of the NLO 
corrections amount to a reduction of —45% compared to 
LO. As far as the scale dependence is concerned a new 
feature appears compared to the three jet rate: The cross 
sections for p = Hj / 4 and p = Hj do not bracket the re- 
sult for p — Hj/2. This is not uncommon for a NLO 



calculation. While the LO cross section is decreasing as 
a function of the renormalization scale — a direct conse- 
quence of the negative value of the beta function — the 
NLO scale dependence may develop a plateau showing 
a flat scale dependence in a restricted region. This is 
in fact what we would consider as an ideal situation. If 
the central scale is close to the extrema two additional 
results obtained by varying the scale by a factor two up 
and down will be both larger or smaller than the cen- 
tral value. This is precisely what happens in case of 
the four jet cross section. As central result we quote in 
Eq. (21 1 the result for p = Hj /2. Assuming that this 



value is already close to the extrema in the restricted 
scale range Hj /4 < p < Hj we set the super-script de- 
scribing the upwards variation to zero. The sub-script 
describing the downwards variation is obtained from 
mm(al mLO (p = # r /4),af eV - NLO (p = #r)). Us- 
ing the explicit results for p — Hj /4 and p — Hj 



„8TeV-NLO 

cf eV - NL % 



(M = 6 T ) 
= #r/4) 



7.91 (0.05) nb, 
4.91(0.15)nb, 



(22) 
(23) 



thus we get the quoted value for the lower scale vari- 
ation band, — 3.24nb. As in the three jet case moving 
from 7 TeV collider energy to 8 TeV increases the cross 
section. For the four jet cross section the effect is with 
about 50% slightly larger than for the three jet rate. We 
also mentioned that the relative numerical uncertainty 
due to the Monte Carlo integration is of the order of 
1%. The larger uncertainty compared to the three jet 
cross section reflects the fact that the integration is much 
more involved. Although not a physical observable it is 
illustrative to study the contributions from different par- 
ton channels. The result is shown in Tab. Q] Note that 
we do not distinguish between quarks and anti-quarks 
in the initial state. For example the qq channel includes 
the contribution from the initial states qq' , qq' ,q?j . Fur- 
thermore we do not distinguish different quark flavours. 
As can be seen the dominant channels are those with 





gg 


qg 


qq 


relative contribution 


37% 


49% 


14% 



Table 1: Relative contribution of the different parton channels 
to pp -> 4 jets at LO. 



one gluon and one quark in the initial state. This is a di- 
rect consequence of the sizeable partonic cross section 
in combination with the largest parton luminosity. 31% 
of the 49% from these channels come from gu — > u + 3g 
and gd — » d + 3g. The largest single process is gg — > 4g 
which contributes 30%. As in the three jet case we show 
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Figure 3: pj distribution in pp — > 4 jets for the leading jet at the LHC 
with a centre-of-mass energy of 8 TeV. 



Figure 4: Rapidity distribution in pp — > 4 jets for the leading jet at the 
LHC with a centre-of-mass energy of 8 TeV. 



in Fig. |3]the pj distribution of the leading jet. The be- 
haviour follows closely what has been observed in the 
three jet case. The scale dependence of the NLO result 
is significantly reduced compared to the LO prediction. 
Again the K-factor is almost constant over a wide pj 
range. We find a negative correction of about 45% — in 
agreement with the findings for the 'inclusive' four jet 
cross section. Only for extreme pj values the K-factor 
changes significantly. For small pj this might again be 
a consequence of soft gluon corrections which would 
require to go beyond fixed order perturbation theory. At 
large pr the uncertainties due to the numerical integra- 
tion become important, so that the results become un- 
reliable. The pj distributions of the remaining jets fol- 
low pretty much the same pattern. The corresponding 
K-factor is even more stable compared to the pj dis- 
tribution of the highest pj jet. The results are shown 



butions of the different jets look remarkable similar. To 
study this in more detail we investigated the ratio 



in Fig. A. 15 Fig. A. 16 and Fig. A. 17 in the appendix. 
In Fig. [4] the rapidity distribution of the leading jet is 
shown. The qualitative findings are similar to the pj 
distribution: Reduction of the NLO cross sections by 
about 45-50%, significant reduction of the scale depen- 
dence and a K-factor which is almost constant over the 
phase space sampled by the distributions. Looking into 
the rapidity distributions of the remaining jets we find 
the same picture. Again we refer to the appendix for the 
corresponding plots. We notice that the rapidity distri- 



Rj = 



do 4 / do 4 
dx\jj dx\{ 



(24) 



The result for Ri is shown in Fig. [5] Two important 
aspects are visible: First of all the ratio is remarkably 
stable with respect to perturbative corrections. Within 
the numerical uncertainties LO and NLO results are in 
perfect agreement. We also show the effects due to scale 
variation. However since in the ratio the leading power 
in a s cancels the scale variation is not necessarily a reli- 
able estimate of the theoretical uncertainty. An alterna- 
tive way to assess the effect of higher order corrections 
would be to compare the ratio expanded in a 4 with the 
naive ratio where we just divide the predictions for the 
numerator and the denominator. However since we use 
a dynamical scale setting procedure this is not possible. 
The second important observation is that to good ap- 
proximation we have R2 ~ 1 consistent with the naive 
expectation. The results for R3 and R4 look very simi- 
lar. Experimentally a measurement of the different ra- 
tios could be used to validate detector efficiencies and 
to further constrain the jet energy scale. 

The stability of the results shown in Fig.[5]is related 
to the fact that in the ratio the almost constant K-factor 
cancels out. In general it might be beneficial to study 



7 



NJet + Sherpa 

pp->4jet at 8 TeV 




-10 1 
Leading jet rapidity 



NJet + Sherpa 

r pp ->• 4 jet at 8 TeV 




400 600 800 

Leading jet pt [GeV] 



Figure 6: Normalized pj distribution of the leading jet in pp — y 4 jets. 

Figure 5: Ratio R2 of the rapidity distribution for the second jet with 
respect to the leading jet. 



normalized distributions since in the ratio many uncer- 
tainties may cancel. This is evidently true for a s but 
should also hold to some extent for uncertainties due to 
the parton distribution functions. In Fig. [6] and Fig. [7] 
we show the normalized distributions for the leading jet. 
Compared to the un-normalized distributions the size of 
the corrections is reduced. In the rapidity distribution 
for example the K-factor becomes close to one. Again 
we stress that the scale uncertainty does not necessar- 
ily provide a reliable estimate of theoretical uncertainty. 
We mentioned in the description of the numerical setup 
that the LO PDFs come with a rather large value for oc s . 
Observing the sizeable NLO corrections — about 40% 
in case of the three jet cross sections and 45% for the 
four jet cross section — one may ask how much of the 
corrections are actually due to the shift in a s . Further- 
more it would be interesting to disentangle the perturba- 
tive corrections of the hard scattering from the change 
of the PDFs when moving from LO to NLO. To do so 
we show in Fig.[8]and Fig.|9]results where we have used 
NLO PDFs together with the respective a s value to eval- 
uate the LO cross sections. We see that in this way the 
size of the NLO corrections is roughly divided by two. 
We also observe that LO and NLO predictions overlap 
taking the scale uncertainty as an uncertainty estimate 
of uncalculated higher orders. From a phenomenologi- 
cal point of view we are lead to the conclusion that using 



the NLO setup in the evaluation of the LO cross sections 
gives a much better approximation to what happens at 
NLO accuracy compared to the default setup with LO 
PDFs. This is a valuable information for the experi- 
mental analysis in cases where the NLO corrections are 
not available or take to long to be evaluated. In gen- 
eral it would be interesting to investigate whether the 
observations we made here in the case of the four jet 
cross section hold true for a wider class of NLO pro- 
cesses. It is conceivable that a similar procedure also 
works for other processes. The reasoning behind the 
large LO CC S value is to account for the NLO corrections 
to the matrix elements which are not considered in the 
LO partonic cross sections. However for this to work 
for a large variety of different processes would require a 
universal K-factor for NLO corrections. Since K-factors 
can be rather different for different processes there is an 
obvious limitation of the procedure. Naively one could 
expect that this becomes in particular true when cross 
sections are considered which are of different order in 
0C S compared to the cross sections which enter PDF fits. 
In such cases using an NLO setting in the evaluation of 
LO cross sections may provide a better approximation 
to the 'full' (NLO) answer. 
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Figure 7: Normalized eta distribution of the leading jet in pp — > 4 jets. 



Figure 8: pj distribution of the leading jet in pp — > 4 jets using NLO 
PDFs in the evaluation of the LO cross section. 



3.4. Comparison of three and four jet production 

As reference value we have calculated also the two 
jet cross section using the same setup as before: 



^8TeV-LO 



,8TeV-NLO 



1234.9(1.2) nb, (25) 
= 1524.9(2.8) nb. (26) 



Combining (not quite consistent in oc s and ignoring soft 



gluon resummation) the results of Eq. (19 1, Eq. (21 



and Eq. (26 1 we estimate the total jet cross section to be 
of the order of 1600 nb. We thus obtain for the 2:3:4 
jet ratios: 1 : 0.05 : 0.005. Only 5% of the multi-jet 
events are three jet events. The four jet topology is fur- 
ther reduced by a factor 1/10. It is interesting to study 
how the ratio between four jet production and three jet 
production behaves as function of the leading jet pj. 
Similar to what has been shown in Ref. [57 1 we analyse 
the ratio 

d04 j do?, 
dpr j dpr 



(27) 



We observe in Fig. 10 that the reduction of the four 
jet cross section with respect to the three jet cross sec- 
tion is mostly due to the low p? region. At large pj 
the fraction of four jet events with respect to three jet 
events raises. In leading order the fraction is about 1/3 
at 800 GeV while in NLO it is close to 1/2. With the 
exception of the low pj region we find that the K-f actor 
is rather constant and takes a value of about 1 .4. Note 
that we have not shown the scale variation in Fig. 10 



Since in the ratio the scale dependence would largely 
cancel the scale dependence will not provide a reliable 
estimate of the uncalculated higher orders. Note that the 
results shown here slightly differ from what has been 
shown in Ref. l57l . We expect the differences to be a 
consequence of the different bin sizes and the different 
R value used in the anti-kt jet algorithm. 



4. Conclusions 

We have presented a study of three and four jet pro- 
duction at the LHC running at a centre-of-mass energy 
of 8 TeV. The virtual corrections were efficiently eval- 
uated using an on-shell unitarity based method imple- 
mented in the publicly available NJet C++ library. We 
have cross checked the virtual corrections for individ- 
ual phase space points as far as possible with exist- 
ing results. We find complete agreement. The com- 
plicated structure of real radiation processes was treated 
within the Catani-Seymour dipole subtraction method. 
We used the Sherpa Monte Carlo program for this. As 
an important cross check of our approach we repro- 
duced results for LHC at 7 TeV published recently by 
the BlackHat collaboration. We find perfect agree- 
ment. We stress that this comparison represents the first 
independent check of the results given in Ref. Il57ll . The 
calculation documented in this article nicely illustrates 
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Figure 10: Comparison of three and four jet production. 



Figure 9: Normalized eta distribution of the leading jet in pp — > 4 jets 
using NLO PDFs in the evaluation of the LO cross section. 



the performance of the NJet library to evaluate QCD 
one-loop corrections. 

For the three and four jet cross sections we find the 
NLO corrections behave similarly to those obtained at 
7 TeV. The NLO corrections reduce the LO cross sec- 
tions by about 40-50% depending on the jet multiplic- 
ity. Compared to the 7 TeV results the cross sections 
are increased by about 50% due to the larger collider 
energy and the related increase in the parton fluxes. We 
have also studied differential distributions. The dynam- 
ical scale setting procedure results in an almost constant 
K-factor leading to very stable results for the normal- 
ized distributions. To pin down the origin of the large 
negative corrections observed in the cross sections and 
in the unnormalized distributions we have analysed the 
impact of the LO PDFs used in the LO predictions. Us- 
ing an NLO setup in the LO calculation shows that the 
large negative corrections are largely due to the change 
in the PDFs when going from LO to NLO with the oc s 
value being of particular importance. 

The corresponding data for all the plots shown in this 



article is tabulated in appendix Appendix B 



We stress that the results presented here are in fixed 
order perturbation theory. Large corrections at low 
pr due to soft gluon emission may require to go be- 
yond fixed order in perturbation theory to achieve a 
good agreement with the data. Significant progress 
has been obtained in the matching of NLO calculations 



with parton shower predictions. For a number of phe- 
nomenologically important processes matched results 
are now available and the automation is pushed forward 
in two different frameworks ((a)MC@NLO 117214741 and 
POWHEG [75] [761). Since the NJet library we devel- 
oped to perform the calculation presented here is pub- 
licly available all required ingredients to include four jet 
production in (a)MC@NLO or POWHEG are available. 
We also note that the NJet library includes the virtual 
corrections required for the calculation of the five jet 
cross section at NLO accuracy. 
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Figure A. 1 1 : pj distribution in pp — > 3 jets for the 2 nd jet at the LHC 
with a centre-of-mass energy of 8 TeV. 



Appendix A. Additional differential distributions 
for three and four jet production 
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Figure A. 12: pj distribution in pp — > 3 jets for the 3 rd jet at the LHC 
with a centre-of-mass energy of 8 TeV. 



In this section we show results for the transverse mo- 
mentum and rapidity distributions of the subleading jets. 
Since they show a behaviour rather similar to the lead- 
ing jet they were not discussed in detail in the main text. 
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Figure A. 13: Rapidity distribution in pp — > 3 jets for the 2 jet at the 
LHC with a centre-of-mass energy of 8 TeV. 
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Figure A. 14: Rapidity distribution in pp — > 3 jets for the 3 rd jet at the 

LHCwithacentre-of-massenergyof8TeV. Figure A.16: p T distribution in pp —» 4 jets for the 3 jet at the LHC 
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Figure A. 15: pr distribution in pp — > 4 jets for the 2 nd jet at the LHC 
with a centre-of-mass energy of 8 TeV. 



Figure A. 17: pj distribution in pp — > 4 jets for the 4 lh jet at the LHC 
with a centre-of-mass energy of 8 TeV. 
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Figure A. 18: Rapidity distribution in pp — > 4 jets for the 2 nd jet at the 
LHC with a centre-of-mass energy of 8 TeV. 
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Figure A. 20: Rapidity distribution in pp — > 4 jets for the 4 th jet at the 
LHC with a centre-of-mass energy of 8 TeV. 



Figure A. 19: Rapidity distribution in pp — > 4 jets for the 3" 1 jet at the 
LHC with a centre-of-mass energy of 8 TeV. 
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PT,i (GeV) 


LO (pb) 


NLO (pb) 










60.0 


— 86.8 


5606.3 (14.6) 


4316.1 (80.8) 


PT.2 (GeV) 


LO (pb) 


NLO (pb) 


86.8 - 


— 117.8 


32747.4 (28.7) 


24813.3 (163.1) 


60.0 


— 86.8 


51490.8 (39.9) 


29576.7 (204.4) 


117.8 


— 153.1 


43942.3 (34.2) 


18293.4 (219.5) 


86.8- 


- 117.8 


43903.6 (30.1) 


25597.6 (174.4) 


153.1 


— 192.6 


25597.3 (22.2) 


13464.8 (116.1) 


117.8 


— 153.1 


18638.5 (17.7) 


10495.8 (97.3) 


192.6 


— 236.3 


11019.8 (13.3) 


6705.3 (49.8) 


153.1 


— 192.6 


7441.0 (9.4) 


4129.0 (38.7) 


236.3 


— 284.3 


4504.2 (6.7) 


2873.1 (21.2) 


192.6 


— 236.3 


2995.0 (5.0) 


1617.8 (19.5) 


284.3 


— 336.5 


1854.8 (3.6) 


1208.5 (15.8) 


236.3 


— 284.3 


1243.1 (3.0) 


660.4 (10.7) 


336.5 


— 392.9 


779.9 (2.1) 


501.9 (5.2) 


284.3 


— 336.5 


522.6 (1.7) 


276.0 (3.6) 


392.9 


— 453.5 


334.4 (1.3) 


220.6 (3.1) 


336.5 


— 392.9 


229.9 (1.2) 


118.9 (2.2) 


453.5 


— 518.4 


145.8 (0.7) 


94.6 (1.5) 


392.9 


— 453.5 


101.6 (0.7) 


54.2 (1.2) 


518.4 


— 587.5 


65.8 (0.5) 


43.2 (0.7) 


453.5 


— 518.4 


46.3 (0.4) 


24.1 (0.7) 


587.5 


— 660.9 


29.2 (0.2) 


19.3 (0.4) 


518.4 


— 587.5 


21.4 (0.2) 


10.9 (0.4) 


660.9 


— 738.5 


13.5 (0.2) 


9.1 (0.3) 


587.5 


— 660.9 


9.7 (0.2) 


4.9 (0.3) 


738.5 


— 820.3 


6.10 (0.09) 


3.9 (0.1) 


660.9 


— 738.5 


4.46 (0.09) 


2.3 (0.1) 


820.3 


— 906.3 


2.84 (0.06) 


2.0 (0.1) 


738.5 


— 820.3 


2.13 (0.06) 


1.2 (0.1) 


906.3 


— 996.6 


1.29 (0.03) 


0.89 (0.06) 


820.3 


— 906.3 


0.98 (0.04) 


0.65 (0.07) 


996.6 - 


- 1091.1 


0.60 (0.02) 


0.49 (0.03) 


906.3 


— 996.6 


0.45 (0.02) 


0.30 (0.03) 


1091.1 


— 1189.8 


0.267 (0.009) 


0.20 (0.01) 


996.6 - 


- 1091.1 


0.22 (0.01) 


0.14 (0.02) 


1189.8 


— 1292.8 


0.123 (0.005) 


0.09 (0.02) 


1091.1 


— 1189.8 


0.099 (0.006) 


0.056 (0.007) 


1292.8 


— 1400.0 


0.054 (0.003) 


0.043 (0.004) 


1189.8 


— 1292.8 


0.043 (0.003) 


0.027 (0.004) 










1292.8 


— 1400.0 


0.020 (0.001) 


0.013 (0.002) 



Table B.2: Table for the p T distribution of the 1 st jet in pp -> 3 jets. 

Table B.3: Table for the p T distribution of the 2 nd jet in pp -> 3 jets. 



Appendix B. Tables of differential distributions for 
three jet production 



For future experimental and theoretical studies it 
might be useful to present the data shown in the figures 
also in numerical form. In the following we give for all 
the plots shown in the article the results obtained in the 
numerical integration. 



PT3 (GeV) 


LO (pb) 


NLO (pb) 


60.0 


— 72.8 


74805.4 (46.5) 


41272.9 (208.2) 


72.8 


— 87.6 


31467.0 (23.7) 


18659.0 (90.8) 


87.6- 


- 104.5 


12270.2 (12.5) 


7475.7 (50.6) 


104.5 


— 123.3 


4819.1 (6.4) 


3008.2 (23.0) 


123.3 


— 144.2 


1928.3 (3.6) 


1244.0 (14.5) 


144.2 


— 167.1 


789.0 (1.9) 


523.2 (6.7) 


167.1 


— 192.0 


326.9 (1.1) 


219.4 (2.7) 


192.0 


— 219.0 


139.4 (0.6) 


95.4 (2.8) 


219.0 


— 248.0 


59.4 (0.3) 


40.6 (0.8) 


248.0 


— 278.9 


26.0 (0.2) 


17.4 (0.7) 


278.9 


— 312.0 


11.48 (0.09) 


8.1 (0.2) 


312.0 


— 347.0 


5.25 (0.06) 


3.7 (0.1) 


347.0 


— 384.0 


2.37 (0.03) 


1.73 (0.09) 


384.0 


— 423.1 


1.10 (0.02) 


0.72 (0.04) 


423.1 


— 464.2 


0.503 (0.010) 


0.34 (0.03) 


464.2 


— 507.3 


0.235 (0.005) 


0.17 (0.01) 


507.3 


— 552.5 


0.107 (0.003) 


0.067 (0.007) 


552.5 


— 599.6 


0.052 (0.002) 


0.041 (0.004) 


599.6 


— 648.8 


0.0244 (0.0009) 


0.015 (0.002) 


648.8 


— 700.0 


0.0113 (0.0005) 


0.007 (0.001) 



Table B.4: Table for the p T distribution of the 3 rd jet in pp -» 3 jets. 
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111 




LO (pb) 


NLO 


(pb) 


-2.8 — 


:-2.5) 


2683.1 (16.1) 


1328.9 


(46.8) 


-2.5 — 


'-2.2) 


3692.9 (15.3) 


1842.1 


(47.8) 


-2.2 — 


'-2.0) 


4669.0 (12.9) 


2477.9 


(37.6) 


-2.0 — 


;-i.7) 


5567.8 (12.0) 


3122.1 


(36.3) 


-1.7 — 


;-i-4) 


6398.3 (11.2) 


3597.6 


(34.4) 


-1.4 — 




7139.7 (10.4) 


4093.4 


(37.6) 


-1.1 — 


;-o.8) 


7739.4 (10.1) 


4613.4 


(90.7) 


-0.8 — 


;-o.6) 


8220.6 (10.0) 


4790.3 


(95.9) 


-0.6 — 


;-o.3) 


8518.6 (9.5) 


5139.7 


(37.9) 


-0.3 — 


(0.0) 


8683.8 (9.5) 


5184.0 


(40.4) 


0.0 — 


0.3 


8696.8 (9.6) 


5208.2 


(51.7) 


0.3 — 


0.6 


8541.3 (9.5) 


5073.1 


(42.3) 


0.6 — 


0.8 


8208.1 (9.7) 


4967.0 


(37.2) 


0.8 — 


1.1 


7730.9 (10.1) 


4571.0 


(34.3) 


1.1 — 


1.4 


7129.7 (10.6) 


4090.0 


(34.6) 


1.4 — 


1.7 


6421.9 (11.3) 


3632.4 


(44.6) 


1.7 — 


2.0 


5580.1 (12.1) 


3113.5 


(45.6) 


2.0 — 


2.2 


4677.1 (13.7) 


2536.4 


(60.3) 


2.2 — 


2.5 


3682.1 (14.8) 


1890.8 


(59.3) 


2.5 — 


2.8 


2670.7 (18.5) 


1299.0 


(42.9) 



Table B.5: Table for the rapidity distribution of the 1 st jet in pp — > 3 
jets. 



T12 




LO (pb) 


NLO 


(pb) 


-2.8 — 


-2.5) 


2815.4 


(13.9) 


1413.7 


(36.2) 


-2.5 — 


-2.2) 


3688.7 


(12.9) 


1921.8 


(35.9) 


-2.2 — 


-2.0) 


4591.1 


(11.6) 


2511.6 


(45.0) 


-2.0 — 


-1.7) 


5473.5 


(11.7) 


3047.2 


(39.3) 


-1.7 — 


-1.4) 


6304.7 


(11.7) 


3578.4 


(48.0) 


-1.4 — 


-1.1) 


7083.9 


(11.9) 


3967.2 


(47.7) 


-1.1 — 


-0.8) 


7727.6 


(11.1) 


4556.2 


(41.9) 


-0.8 — 


-0.6) 


8250.3 


(11.8) 


4901.3 


(46.6) 


-0.6 — 


-0.3) 


8609.4 


(11.5) 


5078.2 


(45.0) 


-0.3 — 


(0.0) 


8803.3 


(11.6) 


5269.1 


(59.3) 


0.0 — 


0.3 


8788.3 


(11.3) 


5284.6 


(46.5) 


0.3 — 


0.6 


8596.2 


(11.6) 


5176.7 


(51.4) 


0.6 — 


0.8 


8234.9 


(11.8) 


4807.0 


(48.9) 


0.8 — 


1.1 


7732.1 


(11.5) 


4462.6 


(45.6) 


1.1 — 


1.4 


7080.3 


(12.9) 


4145.7 


(66.0) 


1.4 — 


1.7 


6319.7 


(12.1) 


3591.5 


(43.4) 


1.7 — 


2.0 


5468.7 


(12.0) 


3067.9 


(37.7) 


2.0 — 


2.2 


4588.0 


(12.5) 


2462.2 


(51.4) 


2.2 — 


2.5 


3680.5 


(12.7) 


1931.2 


(40.3) 


2.5 — 


2.8 


2815.3 


(13.7) 


1396.6 


(46.7) 



Table B.6: Table for the rapidity distribution of the 2 nd jet in pp -> 3 
jets. 



Tl3 




LO (pb) 


NLO 


(pb) 


-2.8 — 


;-2.5) 


3102.7 


(12.3) 


1558.4 


(46.9) 


-2.5 — 


;-2.2) 


3877.2 


(13.3) 


2069.5 


(43.4) 


-2.2 — 


;-2.o) 


4704.9 


(11.2) 


2524.2 


(45.5) 


-2.0 — 


:-i.7) 


5507.7 


(11.5) 


3091.2 


(43.8) 


-1.7 — 


;-i.4) 


6294.6 


(12.4) 


3660.3 


(47.9) 


-1.4 — 




7005.9 


(11.4) 


4030.4 


(42.4) 


-1.1 — 


;-o.8) 


7628.2 


(11.7) 


4385.8 


(51.1) 


-0.8 — 


;-o.6) 


8117.4 


(12.0) 


4766.5 


(40.8) 


-0.6 — 


;-o.3) 


8464.6 


(12.2) 


5028.8 


(42.0) 


-0.3 — 


(0.0) 


8649.2 


(12.9) 


5012.9 


(50.5) 


0.0 — 


0.3 


8652.7 


(12.5) 


5081.1 


(47.0) 


0.3 — 


0.6 


8467.3 


(12.3) 


4897.5 


(54.5) 


0.6 — 


0.8 


8109.7 


(12.0) 


4927.0 


(51.2) 


0.8 — 


1.1 


7626.6 


(12.5) 


4463.0 


(51.9) 


1.1 — 


1.4 


7005.6 


(12.5) 


4112.4 


(68.9) 


1.4 — 


1.7 


6272.6 


(11.4) 


3639.6 


(40.9) 


1.7 — 


2.0 


5487.1 


(11.3) 


3050.4 


(43.4) 


2.0 — 


2.2 


4717.6 


(12.3) 


2577.2 


(55.2) 


2.2 — 


2.5 


3872.4 


(12.1) 


2090.8 


(34.1) 


2.5 — 


2.8 


3087.8 


(12.4) 


1603.9 


(37.7) 



Table B.7: Table for the rapidity distribution of the 3 rd jet in pp -> 3 
jets. 
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PT,i (GeV) 


LO (pb) 


NLO (pb) 


60.0 


92.0 


1013.0 (3.8) 


446.9 (41.5) 


92.0- 


— 129.6 


4204 4 (3 3) 


2335.6 (57.6) 


129.6 


— 172.9 


3775.5 (2.4) 


2303.3 (43.1) 


172.9 


— 221.9 


2626.8 (1.9) 


1452.7 (31.5) 


221.9 


— 276.5 


1498.8 (1.4) 


867.7 (22.6) 


276.5 


— 336.7 


707.1 (0.8) 


417.7 (11.1) 


336.7 


— 402.6 


308.8 (0.4) 


187.0 (5.1) 


402.6 


— 474.1 


131.5 (0.3) 


83.2 (3.1) 


474.1 


— 551.3 


55.4 (0.1) 


32.5 (1.6) 


551.3 


— 634.1 


23.32 (0.08) 


13.8 (0.9) 


634.1 


— 722.6 


9.79 (0.05) 


5.4 (0.5) 


722.6 


— 816.7 


4.09 (0.03) 


2.9 (0.2) 


816.7 


— 916.5 


1.69 (0.02) 


0.9 (0.1) 


916.5- 


- 1021.9 


0.705 (0.009) 


0.40 (0.07) 



Table C.8: Table for the p T distribution of the 1 st jet in pp -> 4 jets. 



Appendix C. Tables of differential distributions for 
four jet production 
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PT.2 (GeV) 


LO (pb) 


NLO (pb) 


60.0 


— 89.8 


3412.4 (4.8) 


2088.7 (58.8) 


89.8- 


- 124.8 


5607.7 (3.1) 


3173.1 (56.4) 


124.8 


— 165.1 


3031.8 (1.9) 


1 1 ( \ A { O A 1 \ 

1670.4 (34.1) 


165.1 


— 210.6 


1342.6 (1.1) 


700.9 (16.3) 


210.6 


— 261.4 


563.4 (0.7) 


304.7 (8.5) 


261.4 


— 317.5 


233.4 (0.4) 


125.9 (4.3) 


317.5 


— 378.8 


97.5 (0.2) 


49.1 (2.3) 


378.8 


— 445.4 


41.5 (0.1) 


21.5 (1.4) 


445.4 


— 517.2 


17.77 (0.08) 


9.0 (0.7) 


517.2 


— 594.2 


7.59 (0.05) 


3.9 (0.3) 


594.2 


— 676.6 


3.27 (0.03) 


1.9 (0.3) 


676.6 


— 764.2 


1.42 (0.02) 


0.7 (0.1) 


764.2 


— 857.0 


0.61 (0.01) 


0.25 (0.09) 


857.0 


— 955.1 


0.264 (0.006) 


0.16 (0.03) 


955.1 - 


- 1058.4 


0.113 (0.004) 


0.06 (0.02) 



Table C.9: Table for the p T distribution of the 2 nd jet in pp -> 4 jets. 



PTj (GeV) LO (pb) NLO (pb) 



Table C. 10: Table for the p T distribution of the 3 rd jet in pp -> 4 jets. 



60.0 


— 76.4 


4670.9 (5.1) 


2710.4 (70.8) 


76.4 


— 95.8 


5196.7 (3.0) 


2941.8 (48.8) 


95.8- 


- 118.0 


2617.1 (1.7) 


1463.5 (25.4) 


118.0 


— 143.2 


1123.2 (0.9) 


619.2 (13.2) 


143.2 


— 171.2 


455.2 (0.5) 


247.8 (6.7) 


171.2 


— 202.2 


180.5 (0.3) 


103.7 (3.0) 


202.2 


— 236.1 


71.1 (0.1) 


37.9 (1.6) 


236.1 


— 272.8 


28.10 (0.08) 


15.9 (0.8) 


272.8 


— 312.5 


11.08 (0.04) 


6.0 (0.4) 


312.5 


— 355.0 


4.46 (0.02) 


2.4 (0.2) 


355.0 


— 400.5 


1.80 (0.01) 


1.0 (0.1) 


400.5 


— 448.9 


0.737 (0.008) 


0.33 (0.05) 


448.9 


— 500.1 


0.281 (0.004) 


0.14 (0.02) 


500.1 


— 554.3 


0.116 (0.002) 


0.064 (0.010) 


554.3 


— 611.4 


0.0434 (0.0008) 


0.017 (0.007) 
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Pta (GeV) 


LU ipo) 




ou.u — 




C\H AC H f C 0\ 

y /43.0 (j.o j 


C A HO A ( OA A \ 

>4/3.4 \6\jA) 


74.2- 


-91.0 


7 1 ~7Q Q (1 (W 

j 1 lo.y [Z.V ) 


lol l.j \Ly .y ) 


91.0 — 


110.2 


984.5 (0.8) 


591.7 (13.6) 


110.2- 


- 131.9 


307.7 (0.4) 


182.8 (5.2) 


131.9 — 


- 156.2 


97.5 (0.2) 


61.2 (2.0) 


156.2- 


- 183.0 


31.49 (0.07) 


19.1 (0.9) 


183.0- 


-212.3 


10.36 (0.03) 


7.2 (0.4) 


212.3- 


-244.1 


3.45 (0.02) 


2.3 (0.2) 


244.1 - 


-278.4 


1.173 (0.008) 


0.79 (0.09) 


278.4- 


-315.2 


0.408 (0.004) 


0.27 (0.03) 


315.2- 


- 354.5 


0.142 (0.002) 


0.12 (0.02) 


354.5 - 


-396.3 


0.0491 (0.0008) 


0.030 (0.006) 


Table C. 1 1 : Table for the p T distribution of the 4 th jet in pp -> 4 jets. 




ill 


LO (pb) 


NLO (pb) 


-2.8- 


-(-2.3) 


490.6 (1.7) 


231.7 (25.6) 


-2.3- 


-(-1.9) 


828.3 (1.8) 


450.5 (24.8) 


-1.9- 


-(-1.4) 


1151.0 (1.6) 


651.9 (22.7) 


-1.4- 


-(-0.9) 


1413.3 (1.8) 


790.3 (24.8) 


-0.9- 


- (-0-5) 


1600.9 (1.7) 


957.7 (21.9) 


-0.5 


-(o.o) 


1698.2 (2.2) 


1005.6 (26.6) 


0.0 


— 0.5 


1698.8 (2.3) 


973.1 (27.1) 


0.5 


— 0.9 


1600.4 (1.6) 


929.5 (23.6) 


0.9 


— 1.4 


1413.9 (1.7) 


811.2 (23.6) 


1.4 


— 1.9 


1152.2 (1.6) 


677.0 (24.1) 


1.9 


— 2.3 


826.0 (1.6) 


440.7 (25.1) 


2.3 


— 2.8 


487.9 (1.6) 


231.2 (32.0) 



Table C.12: Table for the rapidity distribution of the 1 st jet in pp -> 4 
jets. 
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LO (pb) 


NLO (pb) 


-2.8- 


- (-2.3) 


514.8 (1.8) 


262.8 (23.9) 


-2.3- 


-(-1.9) 


819.3 (1.4) 


445.6 (29.8) 


-1.9- 


-(-1.4) 


1126.2 (1.6) 


598.0 (27.7) 


-1.4- 


- (-0.9) 


1399.4 (1.7) 


823.5 (23.4) 


-0.9- 


-(-0.5) 


1602.2 (1.6) 


984.4 (27.8) 


-0.5 


-(o.o) 


1720.2 (2.3) 


1006.7 (23.8) 


0.0 


— 0.5 


1716.3 (2.5) 


972.4 (24.2) 


0.5 


— 0.9 


1607.2 (1.9) 


903.8 (25.9) 


0.9 


— 1.4 


1398.9 (1.8) 


819.4 (26.7) 


1.4 


— 1.9 


1126.3 (1.6) 


628.3 (26.2) 


1.9 


— 2.3 


818.3 (1.5) 


452.8 (22.6) 


2.3 


— 2.8 


512.3 (1.4) 


252.7 (23.4) 



Table C. 1 3: Table for the rapidity distribution of the 2 nd jet in pp -> 4 
jets. 
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LO (pb) 


NLO (pb) 


-2.8 — 


-2.3) 


563.2 (1.3) 


284.2 (21.8) 


-2.3 — 


;-1.9) 


838.6 (1.7) 


472.9 (32.7) 


-1.9 — 


:-i.4) 


1124.0 (1.5) 


644.1 (21.9) 


-1.4 — 


;-o.9) 


1381.6 (1.7) 


787.9 (27.2) 


-0.9 — 


-0.5) 


1582.7 (2.3) 


938.1 (24.7) 


-0.5 — 


(o.o) 


1693.3 (2.7) 


957.9 (26.6) 


0.0 — 


0.5 


1692.3 (1.9) 


1016.1 (27.5) 


0.5 — 


0.9 


1580.8 (1.8) 


928.6 (27.0) 


0.9 — 


1.4 


1382.8 (1.7) 


762.1 (27.4) 


1.4 — 


1.9 


1122.6 (1.6) 


613.0 (25.0) 


1.9 — 


2.3 


838.3 (1.3) 


466.2 (19.5) 


2.3 — 


2.8 


561.3 (1.3) 


279.2 (20.4) 



Table C. 14: Table for the rapidity distribution of the 3 rd jet in pp -> 4 
jets. 
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LO (pb) 


NLO (pb) 


-2.8- 


-(-2.3) 


572.1 (1.5) 


277.2 (17.0) 


-2.3- 


-(-1.9) 


834.6 (1.4) 


451.9 (28.7) 


-1.9- 


-(-1.4) 


1118.4 (1.5) 


652.6 (34.0) 


-1.4- 


- (-0-9) 


1376.9 (1.7) 


803.5 (34.9) 


-0.9- 


-(-0-5) 


1579.3 (2.0) 


917.5 (25.5) 


-0.5 


- (o.o) 


1695.5 (2.5) 


989.0 (27.9) 


0.0 


— 0.5 


1695.1 (2.0) 


1002.0 (30.4) 


0.5 


— 0.9 


1585.7 (2.4) 


916.9 (28.1) 


0.9 


— 1.4 


1379.1 (1.8) 


789.6 (26.3) 


1.4 


— 1.9 


1115.4 (1.5) 


638.2 (24.8) 


1.9 


— 2.3 


837.7 (1.5) 


428.9 (20.0) 


2.3 


— 2.8 


571.7 (1.1) 


283.1 (18.9) 



Table C. 15: Table for the rapidity distribution of the 4 th jet in pp -> 4 
jets. 
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